Coherent excitation-energy transfer and quantum entanglement in a dimer 
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We study coherent energy transfer of a single excitation and quantum entanglement in a dimer, which consists 
of a donor and an acceptor modeled by two two-level systems. Between the donor and the acceptor, there exists a 
dipole-dipole interaction, which provides the physical mechanism for coherent energy transfer and entanglement 
generation. The donor and the acceptor couple to two independent heat baths with diagonal couplings that do 
not dissipate the energy of the non-coupling dimer. Special attention is paid to the effect on single-excitation 
energy transfer and entanglement generation of the energy detuning between the donor and the acceptor and the 
temperatures of the two heat baths. It is found that, the probability for single-excitation energy transfer largely 
depends on the energy detuning in the low temperature limit. Concretely, the positive and negative energy 
detunings can increase and decrease the probability at steady state, respectively. In the high temperature limit, 
however, the effect of the energy detuning on the probability is negligibly small. We also find that the probability 
is negligibly dependent on the bath temperature difference of the two heat baths. In addition, it is found that 
quantum entanglement can be generated in the process of coherent energy transfer. As the bath temperature 
increases, the generated steady-state entanglement decreases. For a given bath temperature, the steady-state 
entanglement decreases with the increase of the absolute value of the energy detuning. 
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PACS numbers: 03.65.Yz, 71.35.-y, 03.67.Mn 



I. INTRODUCTION 

Coherent excitation energy transfer is an important step of 
photosynthesis 1 1], in which photosynthetic pigments capture 
the solar light to create electronic excitations and then trans- 
fer the excitation energy to a reaction center fZ-'T]. Usually, 
the transfer of a single excitation from the pigment where the 
electronic excitation is created to the reaction center is a very 
complicated physical process, since the practical transfer pro- 
cess takes place on a complicated network of pigments. How- 
ever, the basic physical mechanism can be revealed in such 
a light-harvesting complex by studying a basic part: a dimer 
system which consists of a donor and an acceptor modeled by 
two two-level systems. 

On a complicated network of pigments, there generally ex- 
ist two kinds of interactions. On one hand, between any two 
pigments there exists a dipole-dipole interaction, which re- 
sults in excitation energy transfer. On the other hand, the 
pigments interact inevitably with their surrounding environ- 
ments such as the nuclear degrees of freedom and the proteins. 
Corresponding to different cases for the scale of the two kind 
of interactions, different approaches have been proposed to 
study the single-excitation energy transfer Concretely, when 
the dipole-dipole interactions between any two pigments are 
much weaker than the interactions of the pigments with their 
environments, the energy transfer process can be well char- 
acterized by the Forster theory |8], in which the evolution of 
the network is calculated perturbatively up to the second order 
in the dipole-dipole interactions between the pigments; When 
the interactions of the pigments with their environments are 
much weaker than the dipole-dipole interactions between any 
two pigments, various approaches based on the quantum mas- 
ter equation have been proposed (e.g., Refs. |9-23]), in which 
the evolution of the network is calculated perturbatively up to 



the second order in the interactions between the pigments and 
their environments. 

With the above considerations, in this article we study 
single-excitation energy transfer in a dimer, which consists of 
a donor and an acceptor modeled by two two-level systems. 
Obviously, when the donor and the acceptor are decoupled, it 
is impossible to realize energy transfer between them. There- 
fore, the simplest way to realize energy transfer is to turn on 
a non-trivial interaction (for example, the dipole-dipole inter- 
action) between the donor and the acceptor Then a single 
excitation can coherently oscillate between the donor and the 
acceptor. However, in this case, there is no steady-state energy 
transfer, namely the transferred energy can not approach to a 
stationary value. In the presence of environments, the donor 
and the acceptor will inevitably couple with environments. In 
general, the coupling form between the donor (acceptor) and 
its environment is diagonal in the representation of the free 
Hamiltonian of the donor (acceptor). Physically, due to this 
type of coupling, although the excitation energy will not de- 
cay into the environments, it will induce a steady-state energy 
transfer between the donor and the acceptor Since in practi- 
cal cases both the characteristic frequency and the heat bath 
temperatures of the donor and the acceptor may be different 
due to different chemical structures, we study in detail how 
the characteristic frequencies and the heat bath temperatures 
of the donor and acceptor affect the efficiency of the excita- 
tion energy transfer This is one point of the motivations of 
our present work. 

In the presence of the interactions between the pigments 
for transferring energy, a naturally arising question is how 
about the quantum entanglement among the pigments which 
are involved in the energy transfer process. Because quan- 
tum entanglement is at the heart of the foundation of quantum 
mechanics |2J, 12511 and quantum information science (e.g.. 



Refs. 11261 12711 ). it is interesting to know how is the dynam- 
ics of the created quantum entanglement in the dimer sys- 
tem during the process of single-excitation energy transfer 
This is the other point of the motivations of our present in- 
vestigations. In fact, recently people have become aware of 
quantum entanglement in some chemical and biologic sys- 
tems (e.g., Refs. 1 17, 28-32]) such as photosynthetic light- 
harvesting complexes |17, 31, 32]. 

This article is organized as follows: In Sec. [Ill we present 
the physical model and the Hamiltonian for studying the 
single-excitation energy transfer. A dimer consists of a donor 
and an acceptor, which are immersed in two independent heat 
baths. Between the donor and the acceptor, there exists a 
dipole-dipole interaction, which provides the physical mecha- 
nism for coherent excitation energy transfer and entanglement 
generation. In Sec. Hill we derive a quantum master equation 
to describe the evolution of the dimer. Based on the quantum 
master equation we obtain optical Bloch equations and their 
solutions. In Sec. lIVI we study single-excitation energy trans- 
fer from the donor to the acceptor. The effect on the transfer 
probability of the energy detuning and the bath temperatures 
are studied carefully. In Sec. |V] we study the quantum en- 
tanglement between the donor and the acceptor by calculating 
the concurrence. We conclude this work with some remarks in 
Sec. [yi] Finally, we give an appendix for derivation of quan- 
tum master equation (|7]l. 



II. PHYSICAL MODEL AND HAMILTONIAN 

As illustrated in Fig. [IJa), the physical system under our 
consideration is a dimer, which consists of a donor and an 
acceptor modeled by two two-level systems (TLSs), TLSl 
(donor) and TLS2 (acceptor), with respective energy separa- 
tions oji and W2- The donor and the acceptor are immersed 
in two independent heat baths of temperatures Ti and T2, re- 
spectively. Between the donor and the acceptor there exists 
a dipole-dipole interaction of strength ^. The Hamiltonian of 
the total system, including the two coupled TLSs and their 
heat baths, is composed of three parts. 



H — H-T 



. + Hn+ H, 



(1) 



where //tlss is the Hamiltonian (with h - 1) of the two cou- 
pled TLSs, 






Concretely, the first two terms in Eq. ^ are free Hamil- 
tonians of the two TLSs, which are described by the usual 
Pauli operators cr,+ = (crp"'' - {cr_^ + /o"v)/2 = \e)ii {g\ and 
cTi - \e)ii {e\ - \g)ii {g\, where \g)i and |e)/ are, respectively, the 
ground and excited states of the /th (/ = 1,2) TLS, namely 
TLS/. The last term in Eq. (|2]i depicts the dipole-dipole inter- 
action of strength ^ between the two TLSs. This dipole-dipole 
interaction provides the physical mechanism for excitation en- 
ergy transfer and entanglement generation between the two 
TLSs. 
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FIG. 1: (Color online) (a) Schematic of the physical system. A 
donor and an acceptor are immersed in two independent heat baths 
of temperatures Ti and Ti, respectively. A dipole-dipole interac- 
tion of strength ^ exists between the donor and the acceptor, which 
are described by two two-level systems with resonant frequencies aii 
and u>2, respectively, (b) The energy levels of the bare states \r]„) 
(n = 1,2, 3,4) of the donor and the acceptor when they are decou- 
pling, (c) The energy levels of the eigenstates \A„) (n = 1, 2, 3, 4) of 
the coupled donor and acceptor. The corresponding eigen-energies 
are denoted by E„. The parameters r23 and r32 are, respectively, the 
bath induced transition rates from states |/1t) to \At,) and from states 

l^3> to 1^). 



The Hilbert space of the donor and the acceptor is of four 
dimension with the four basis states Irji) - \ee), I772) - \eg), 
Vli) - \g^)^ and I/74) = \gg), as shown in Fig. [TJb). In the 
presence of the dipole-dipole interaction, a stationary single- 
excitation state should be delocalized and composed of a com- 
bination of the single-excitation in the two TLSs. According 
to Hamiltonian (|2]i, we can obtain the following four eigen- 
states 



M2) 
1^3) 
1^4) 



\ee), 

cos.{9l2)\eg) + sm(6l2)\ge) , 

-sin(0/2)|e^) + cos(0/2)|ge). 



(3) 



(2) and the corresponding eigenenergies/ii = -E4 - (co\ +0)2) 12 
and £2 = -E-i - ^JAaJ^/4 + ^^, as shown in Fig.ITJc), by solv- 
ing the eigen-equation //tlss l^n) - E„\An) (n - 1,2,3,4). 
Here we introduce the energy detuning Aw = oji - a;2 and the 
mixing angle defined by tanO = 2^/ Aw. Note that here 
the mixing angle < 6 < n. Therefore, when Aw > 0, 
namely wi > W2, we have = arctan(2^/Aa;); however, when 
Aw < 0, that is wi < W2, we have - arctan(2f/Aw) + n. 

As pointed out by Caldeira and Leggett pS], when the cou- 
plings of a system with its environment are weak, it is uni- 
versal to model the environment of the system as a harmonic 
oscillator heat bath. In this work, we suppose that the cou- 



plings of the TLSs with their environments are weak, then it 
is reasonable to model the environments as two harmonic os- 
cillator heat baths with the Hamiltonian 



picture reads. 
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(4) 
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Here H'g' and Hg 



m 



are respectively the Hamiltonians of the 
heat baths for the TLSl and TLS2, 



(a) 



H' 



,t /A-N 



= 2 



OJaj^ -^ }•> 



Hf ^Yj'^bkblbk, (5) 



where a. {b[) and Oj {bk) are, respectively, the creation and an- 
nihilation operators of the jth (^th) harmonic oscillator with 
frequency ojaj (cobk) of the heat bath for TLSl (TLS2). In 
practical systems of excitation energy transfer, the environ- 
ment is composed of the nuclear degrees of freedom of the 
molecules. 

The interaction Hamiltonian of the TLSs with their heat 
baths reads (e.g., Refs. I.9.-.17.1) 

H! = o-[cr7 2] 8iM] + Clj) + 0-2 0"2 /_j 82k(bl + bk). (6) 



In this case, there is no energy exchange between the TLSs 
and their heat baths. This type of diagonal coupling has been 
used to describe the dephasing of quantum systems i34ll . For 
simplicity, but without loss of generality, in the following we 
assume the coupling strengthes gij and g2k are real numbers. 



III. QUANTUM MASTER EQUATION AND OPTICAL 
BLOCH EQUATIONS 

Generally speaking, there are two kinds of different ap- 
proaches to study photonsynthetic excitation energy transfer 
One is based on the Forster theory [8], which is valid when 
the electronic couplings between pigments are smaller than 
the couplings between electrons and environments. ^he other 
is usually based on quantum master equations |9-23] in var- 
ious forms, which are valid when the electron-environment 
couplings are weaker than electronic couplings between pig- 
ments. In this work, we shall consider the latter case where the 
coupling (with strength ^) between the two TLSs is stronger 
than the couplings (relating to y) between the TLSs and their 
local environments (in our following considerations we take 
^/y = 5). We will derive a quantum master equation by 
truncating the evolution up to the second order in the TLS- 
environment couplings. On the other hand, we derive the mas- 
ter equation in the eiegen-representation of the two cou pled 
TLSs so we may safely make the secular approximation 133] 
by neglecting the high-frequency oscillating terms. This ap- 
proximation is also equivalent to rotating wave approximation 
in quantum optical systems. The detailed derivation of the 
quantum master equation will be presented in the appendix. 

In the eigen-representation of Hamiltonian (|2]) of the two 
coupled TLSs, the quantum master equation in Schrodinger 



+ 2^ '^n{'^0-n„PsCr„n-Cr„„ps - PsCTnn) 
n=l,2,3 

+r32(2cr23PsO"32 - 0"33PS - PsCT^i) 

+r23(2o-32PsCr23 - O-22PS -P5 0"22) 

+2X12(0-1 IPs 0-22 -I- cr22pscrii) 

+2X13(0-1 IPs 0-33 + CTiipsO-u) 

+2X23(o-33PsCr22 + cr22PsO-33). 



(7) 



In Eq. ([7|i, ps is the reduced density matrix of the two TLSs. 
The transition operators cr„,„ (n, m = 1, 2, 3, and 4) are defined 
as o-,„„ = \A„){A,„\, where the states \A„) have been defined in 
Eq. (O. Meanwhile, we introduce the effective rates as fol- 
lows; 



Hi 


= X\ +X2^ 


n2 


= co^\9l2)Xx + ^m\0l2)x2. 


n3 


= sin^0/2ki + co^\ei2)x2. 


r32 


- - sin^0[yini(e) -i-'X2n2(e)] 


r23 


= -sin^0[yi(ni(e)-H l) + 72( 


Xn 


= cosHemxx + ^in\ei2)x2. 


Xu 


= ^m\ei2)xy + co^\ei2)x2. 


Xii 


= -sin^ecri +X2), 



(8) 



where xi = lim^^0'S;(w)[2n/(w) -Hi], with S i(u) = 
7Tgi(aj)g^(cj) and y; = ngi(E)g^{E) for Z - 1,2. Here gi(aj) 
and g2(oL>) are respectively the densities of state for the two in- 
dependent heat baths surrounding the donor and the acceptor 
The parameter e = £2 - £3 is the energy separation between 
the two eigenstates I/I2) and lA^,}. And 



ni(a>) = 



1 



exp(a> /Ti) - 1 



(9) 



is the thermal average excitation numbers of the heat baths 
of TLS/. Hereafter we set the Boltzmann constant ^^ = 1. 
We consider a special case of the ohmic spectrum densities 
S i{oj) = T]\a) and 82(0^) - ri2CL>, and then we obtain ;t'i = 
2r]iTi and;t'2 - 27/2^2. 

From quantum master equation (|2|, we can see that there 
exist both dissipation and dephasing processes in the eigen- 
representation of the Hamiltonian (|2]). The first line in Eq. (|7| 
describes the unitary evolution of the system under the Hamil- 
tonian (|2|. The second line in Eq. (|7|i describes the dephasing 
of the states \Ai), I/I2), and I/I3). The third and fourth lines de- 
scribe, respectively, the exciting process from I/I3) to I/I2) and 
the decay process from I/I2) to I/I3), as illustrated in Fig.[Tlb). 
Moreover, there exist three cross dephasing processes in the 
last three lines in Eq. (|7]), these terms can decrease the coher- 
ence between two levels, which can be seen from the follow- 
ing optical Bloch equations ( fTOb . 



According to quantum master equation ([7]l, we can de- 
rive optical Bloch equations for the elements {cr„j„{t)) - 

Trs [pXf)0"mn], 

i&nit)) = <(r44(0) = 0, 

<0-22 (0) = - <d-33 (0) = 2r32 (0-33 (0) - 2r23 <0-22 (0) , 

<o-32 (0) = [-is - (n2 + n3 + r23 + r32 - 2x23)] (0-32 (0) • 

(10) 

Here we present only the equations of motion for the elements 
which will be used below. In fact, the equations of motion for 
all of the elements in the density matrix ps can be obtained ac- 
cording to quantum master equation (|7]). Clearly, from optical 
Bloch equations ( fTOl ) we can see that the diagonal elements 
decouple with the off-diagonal elements. It is straightforward 
to get the transient solutions of optical Bloch equations (flOb . 



{<Tn (0) 
{a-12 (0) 



<O-33(0> = 



{(r,2{t)) 



{(Turn, <(r44(0> = <(^44(0)>, 

«cr22(0)) + <cr33(0)))r32 



r23 + r- 



-2(r23+r32)/ 



32 
^ «Cr22(0))r23-<tr33(0))r32) , 

r23 + r32 

«cr22(0)) + <tr33(0)))r23 

r23 + r32 

^ «Cr33 (0)) r32 - <tr22 (0)) r23) _^-2(r„+r„)r 

r23 + r32 

<cr32 (0)> e-(r23+r32+cos= en.y^-iet^ 



(11) 



Here we have used the relation 112 + 03 - 2X23 
The steady-state solutions of Eq. ( fTTT ) read 



cos^ en 



(o"li(oo)) 
(0"22 (00)) 

(0-33 (0°)) 
(0"32 (00)) 



<(r,i(0)), <(r44 (cx))) = ((^44 (0)) , 

«tr22(0)) + <^33(0)))r32 

r23 + r32 

«(r22(0)) + <cr33(0)))r23 



r23 + r32 



0. 



(12) 



The steady-state solutions for other off-diagonal elements of 
the density matrix are zero. Therefore, we can see that the 
steady state of the two TLSs is a completely mixed one. 



IV. PROBABILITY FOR SINGLE-EXCITATION ENERGY 
TRANSFER 

In order to study the probability for single-excitation energy 
transfer from the TLSl (donor) to the TLS2 (acceptor), we 
assume that the TLSl initially possesses a single excitation 
and the TLS2 is in its ground state, which means the initial 
state of the two TLSs is 



1^(0)), 



cos (0/2) 1^2) -sin (0/2) 1^3). (13) 



Since the couphngs between the TLSs and their heat baths are 
diagonal, there is no energy exchange between the TLSs and 
their heat baths, and the probability for finding the TLS2 in 



its excited state is right that of the single excitation energy 
transfer, 

P(t) = Trzlpicr^cr^] 

= {crn(t)) + sm\e/2){cr22(t)) + cos\e/2){cr33(t)) 
+ sin0Re[<cr23(O)], (14) 

where p2 = Tri [ps ] is the reduced density matrix of the TLS2. 

A. Transient state case 

According to Eq. ( fTTI ). the probability given in Eq. (fT4b can 
be expressed as follows: 



P(t) = 



r32sin^(e/2) + r23Cos^(g/2) 

r23 + r32 

+ ,0, /32sin^(e/2)-r23Cos^(e/2) _^_,,r„.r„v 



r23 + r32 



Um^0cos{st)e-(^'-'^^''-^"''"'''^>. 



(15) 



Now, we obtain the probability for single-excitation energy 
transfer from the TLSl to TLS2. This probability ( fTsT i is a 
complicated function of the variables of the two TLSs and 
their heat baths, such as the energy separations oji and W2, the 
strength ^ of the dipole-dipole interaction, and the tempera- 
tures Ti and T2 of the heat baths. To see clearly the effect on 
probability (fTST i of the bath temperatures and the energy sep- 
arations of the TLSs, we introduce the following variables: 
mean temperature T„ = (Ti + T2)/2, mean energy separation 
w,,, = (wi -I- W2)/2, temperature difference AT - T\ - T2, 
and energy detuning Aoj - coi - W2. And Aw > and 
Acij < mean the positive and negative detunings, respec- 
tively. For simplicity, in the following considerations we as- 
sume j\ -j\- J. 

In the following we consider three special cases: (1) The 
resonant case, in which the two TLSs have the same energy 
separations, i.e., (X)\ - (x)2 - cOm, that is Aw = 0. Now the 
mixing angle 6 = n/2 and the energy separation s - 2^. From 
Eq. (fT5] l we obtain 



/^res(0=^-^COS(2^0e4A'(2Or', 



where we introduce the parameter 

A?(2^) = «i(2^) + «2(2^)+l. 



(16) 



(17) 



The subscript "res" stands for resonant case. Equation ( fTSl ) 
means that the probability Pres increases from an initial value 
to a steady-state value 1/2 as the time t increases. How- 
ever, the increase of the probability is exponential modulated 
by a cosine function rather than monotone. In the short time 
limit it may experience small oscillation. The exponential rate 
N{2^)yl2 is a function of the parameters ^, y, T\, and T2. Ob- 
viously, the parameter A^(2^) increases with the increase of the 
temperatures of the heat baths. In the low temperature limit, 
i.e., Til{2^) X and r2/(2^) * 0, we have ni(2^) x and 
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FIG. 2: (Color online) The probability Pres given in Eq. (O is plotted 
vs the scaled time yt for different bath temperatures T„Jy = 0.1 
(solid red line), 10 (dash dotted blue line), and 100 (dashed black 
line) in the resonant case Aco/y = 0. Other parameters are set as 
y = 1, ^/y = 5, and AT/y = 0. 



FIG. 3: (Color online) The probability Phti given in Eq. J19t vs the 
scaled time yt for different mixing angle 6 = 0.6n (solid red line), 
O.Stt (dash dotted blue line), and 0.9ff (dashed black line) at the high 
temperature limit T„Jy = 100. Other parameters are set as y = 1, 
^7 = 5,xdy=Xily = O.Oir„, and ATly = 0. 



n2(2^) ~ 0, then N{2^) * 1. On the contrary, in the high tem- 
perature limit, i.e., Ti/(2^) » 1 and T2/(2^) » 1, we have 
ni(2^) * ri/(2^) and n2(20 ^ T2l{2^\ then 






(18) 



The above equation means that in the high temperature limit, 
the rate N{2^) is proportional to the mean temperature T„, and 
does not depend on the temperature difference IsJ . In Fig.|2l 
we plot the probability Pres vs the scaled time yt for different 
bath temperatures T^, here we assume that T\ - T2 - I'm- 
From Fig. |2] we can see that in the low temperature limit 
the probability increases with an initial oscillation. With the 
increase of the bath temperatures, the oscillation disappears 
gradually. 

(2) The high temperature limit case, i.e., T\,T2 » e. In 
this case, ni(e),n2(e) » 1, then we can make the approxima- 
tions «i(e) ~ «i(e) + 1 and n2(e) ~ "2(6) + 1, which lead to 
r23 ~ r32. Therefore from Eq. (fTSl l we can obtain the time 
dependent probability 



^hti it) 



i-icos2 0e-^'"'«'V(^'^' 
2 2 



1 . 
-— sin" ycos 

2 



(fiOe"*-^"""' "^^^ "'"^ eN{s)y)t^ j-jg-j 



where we introduce the parameter A^(e) = hiie) + h2(s) + 1 
and the subscript "htl" stands for the high temperature limit. 
Obviously, the above probability Phti increases from an initial 
value to a steady-state value 1/2. And the increase of Fhti is 
not simply exponential. In Fig. [3] we plot the probability Pm 
vs the scaled time -yt and the mixing angle in the high tem- 
perature limit. Since the probability ( fT9l ) is a function of sin^ 6 
and cos^ 0, therefore in Fig. [3] we only need to plot the prob- 
ability in Eq. ( fT9] l for the negative detuning cases. Figure [3] 
shows that in the long time limit the probability reaches 1 /2 ir- 
respective of the 9. Note that here the mixing angle Q < 6 < n. 



The cases of < < 7r/2 and 7r/2 < 6 < n mean the energy 
detuning Aw > and Aw < 0, respectively. And the angle 
- nl2 corresponds to the resonant case. Here we choose 
Q.ln < < 0.971, which corresponds to 6.2 > Aw/^ > -6.2. 

(3) The low temperature limit case, i.e., T\,T2 ~ 0. Now 
we can make the approximations ni(e) ~ and «2(e) ~ 0, 
which lead to r32 ~ and r23 ~ sin^ y/2. Then we obtain 
the probability 



Piti(0 ~ cos2(6»/2)(l - COS 0e- ''"''''") 



— sin2 6'cos(ef)e-5™-er' 
2 



(20) 



where the subscript "Itl" means the low temperature limit. In 
this case, the probability increases for an initial value to a 
steady-state value cos^(6'/2). In Fig.JH we plot the probabil- 
ity Piti vs the time t for different mixing angles 6 in the low 
temperature limit. Figure |4] shows that the probability Piti in- 
creases from to a steady state value with the increase of the 
time t. In the short time, the probability experiences small os- 
cillation. The steady state value decreases with the increase of 
the 0. Actually, the obtained results are very reasonable from 
the viewpoint of energy conservation. For the case of < 7r/2, 
the energy detuning Aw > 0, we have u)\ > cl>2, then the en- 
ergy emitted by TLSl can excite more than one TLS2 into 
their excited state; For the case of > 7t/2, we have Aw < 0, 
we have wi < W2, then the energy emitted by TLSl can only 
excite less than one TLS2 into the excited state. Therefore, it 
is understandable that the steady-state value of probability in 
low temperature increases as the parameter decreases. 
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FIG. 4: (Color online) The probability Pui(t) given in Eq. ilOl vs 
the scaled time yt for different mixing angle = OAn (solid red 
line), OAn (dashed brown line), 0.6n (dash dotted blue line), and 
0.9n (solid black line) at the low temperature limit T„,/y = 1. Other 
parameters are set as y = 1, f /y = 5, and AT/y = 0. 



FIG. 5: (Color online) The steady state probability P,,, vs the bath 
temperature T,„ for different mixing angle 6* = 0.1;r (dashed red line), 
0.5 (solid black line), and 0.9n (dash dotted blue line). Other param- 
eters are set as y = 1, ^/y = 5, and AT/y = 0. 



B. Steady state case 



At steady state, the probability ( fTSl l becomes 



1 / cose* 
"' 2 \ N(e) 



(21) 



where the subscript "ss" stands for steady state and A^(e) - 
hi(e) + h2(s) + 1. This steady-state probability is a very in- 
teresting result since it depends on the mixing angle and the 
bath temperatures Ti and T2 independently. It depends on the 
mixing angle 6 and the bath temperatures T[ and T2 by cos 
and 1 /N(s), respectively. 

We first consider several special cases at steady state: (1) 
The resonant case, i.e., Aoj - 0. In this case, cos9 = 0, 
then Pss - 1/2. In the resonant case, the steady-state prob- 
ability P55 for single-excitation energy transfer is indepen- 
dence of the temperatures of the two heat baths. This result 
can also be understood from the following viewpoints: When 
sin(0/2) = cos(0/2) - 1/ V2, the eigenstates I/I2) and I/I3) be- 
come \A2} = (leg) + Ige))/ V2 and 1^3) = Heg) + \ge))/ V2. 
Therefore for any statistical mixture pss = picrii + PiO-^i of 
the two eigenstates I/I2) and \At,), the probability for finding 
the two TLSs in state \ge) is 1/2, where p2 + p^, - I is the 
normalization condition. (2) The high temperature limit, i.e., 
Ti,T2 » e. In this case, «i(e) » 1 and naC^) ^ 1. there- 
fore N{e) » 1, which leads to Pss ~ 1/2. In fact, in the 
high temperature limit, the steady state of the TLSs should be 
Ps ~ (0-22 -I- cr33)/2, therefore according to Eq. (O we know 
that the probability for finding the two TLSs in state \ge) is 
1/2. (3) The low temperature limit, i.e., T\,T2 <^ s. In this 
case, ni(e) ~ and n2{s) » 0, then N{s) » 1, which means 
Pss = cos^{6/2). In Fig.|5] we plot the steady-state probability 
Pss in Eq. ( |2TI ) vs the bath temperatures T,„. Figure |5] shows 
that, for the positive detuning case, i.e., Q < < n/l, the 
steady state probability Pss decreases from 1 to 1/2, but for 
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FIG. 6: (Color online) The steady state probability Pj, vs the mix- 
ing angle 9 for different bath temperature T„,/y = 0.1 (dashed red 
line), 10 (dash dotted blue line), and 100 (solid black line). Other 
parameters are set as y = 1, ^/y = 5, and AT/y = 0. 



the negative detuning case, i.e., n/2 < < tt, the Pss increases 
from to 1/2. For the resonant case, the Pss is 1/2 iiTespec- 
tively of the bath temperature Ti - T2 - T,„. In Fig. |6] we 
plot the steady-state probability Pss in Eq. ( 1211 1 vs the mixing 
angle 0. Figure |6] shows that, in the high temperature case, the 
Pss becomes approximately a fixed value 1 /2 irrespective of 
the 0. But in the low temperature case, the steady state prob- 
ability Pss decreases with the increase of 0. These results are 
consistent with the above analysis. Therefore, in the low tem- 
perature limit, we can improve the steady-state probability Pss 
via increasing the Aw. 

In the above discussions of the steady-state probability, we 
have assumed the bath temperature difference Ar is zero. Ac- 
tually, we also study the dependence of the steady-state prob- 
ability on the bath temperature difference Ar in both the low 
and the high temperature limits. We found that the depen- 



dence of the probability on AT" is negligibly small with the 
current parameters. This result is well understood from the 
following viewpoint: in the low temperature limit, we have 
Ti, T2 <s s, therefore ni(e) ~ and n2(e) ~ 0, N{s) ^ 1, then 
Pss = cos^(0/2), which does not depends on the bath tempera- 
ture difference AT"; On the other hand, in the high temperature 
limit, T\,T2 » e, therefore ni(e) » 1 and h2{s) » 1, then 



bare state representation and the evolution of the system is ex- 
pressed in the eigenstate representation. Therefore we need 
to obtain the transformation between the two representations. 
The density matrix elements in the eigenstate and bare state 
representations are expressed by (cr,j(f)) and {Tij{t)), respec- 
tively. Making using of Eq. Q, we can obtain the relations 
for diagonal density matrix elements 



1 / ecosfl 

Z \ Zi HI 



(22) 



which is independent of the bath temperature difference AT". 



QUANTUM ENTANGLEMENT BETWEEN THE 
DONOR AND ACCEPTOR 



<0-ii(f)> = <Tii(0), (0-44(0) = <T44(0> , 

(cr22(0) = C0S2 (0/2) <T22(0> + Sin2 {612) <T33(0> 

+ -sin0((T23(O> + <T32(f)», 

<cr33(0) = sin^ (612) <T22(0> + cos^ (0/2) <T33(0) 



--sin0((T23(f)) + <T32(f))), 



(28) 



In this section, we study the quantum entanglement be- 
tween the donor and the acceptor with concurrence, which 
will be defined below. For a 2 x 2 quantum system (two TLSs) 
with density matrix p expressed in the bare state representa- 
tion, its concurrence is defined as IBol 



(23) 



dp) = max{0, V^ - V«2 - V«3 - V«4), 



where s,- (/ = 1,2,3,4) are the eigenvalues {s\ being the largest 
one) of the matrix pp, where the operator p is define as 



p = (cr. ® cr\)p*{(r. ® crl) 



(24) 



with p* being the complex conjugate of p. Note that here cr'. 
is the usual Pauli matrix pointing the y axis. For the 2x2 
quantum system, the concurrences C -Q and C - I mean the 
density matrix p is an unentangled and maximally entangled 
states, respectively. Specially, for the "X"-class state with the 
density matrix 



P = 



( pn pi4 

P22 P23 

P32 P33 

P41 P44 ; 



(25) 



expressed in the bare state representation, the concurrence 
is 03711 

C{p) = max{0,2(|p23l - VPllP44),2(|pi4| - VP22P33))- (26) 

Now, for the present system, its density matrix p can be ex- 
pressed as the following form in the bare state representation. 



and the following off-diagonal element which will be useful 
below. 



{(Tliit)) = :^^ine{{T^^{t)) - {T22{t))) 

+ cos^ (6/2) <T23(0> - sin^ (6/2) <T32(0) -(29) 



Correspondingly, we can obtain the inverse transform 



<T22(0> = cos^ (6/2) <cr22(0) + sin2 (6/2) <cr33(f)) 

-- sin 6 ({0-23(1)) + {0-32(1))) , 
<T33(0) = sm^(6/2){cr22(t)} + cosH6/2){cr33(t)) 

+ - sin 6 ({o-23(t)) + {0-32(1))) , 
<T23(0) = - sin^ (6/2) {(r32(t)} + cos^ (6/2) {cr23(t)} 

+ ^sin6({cr22(t))-{cr33(t))). (30) 



Also here we only express explicitly the elements which will 
be used below. 

In order to calculate the concurrence of the system, we 
need to know its density matrix in the bare representation for 
a given initial state. Fortunately, the evolution relation from 
(Ty(0)) to {Tjj(t)} can be obtained through the following pro- 
cess 



({tu) <T2I> <T3i) <T4i)^ 

<Tl2) <T22) <T32) <T42> 

<Tl3) <T23) <T33) {T43} 

<Tl4> <T24) <T34) <T44> j 



(27) 



where the density matrix elements are defined as {t,j) - 
Tr[T,jp] = Tr[|77,)(?7j|p] = {'Hjlpl^i) with the transition oper- 
ator Tij = \rii}{T]j\. Since the concurrence is defined in the 



<T,/0)) ^ {cTij(0)) -^ {cTijd)) -. <r,/0). 



(31) 



Concretely, the transformation relations (ry(0)) — > (cr,-y(0)) 
and {o-ij(t)} — » {Tij(t)} are determined by Eqs. ( |28] ), ( |29l ) 
and ( l30l l, and the evolution relation (cr,y(0)) — > {a-ij(t)) is de- 
termined by Eq. (HUi. In terms of Eqs. ^, ^, ^, ( |30] |. 
and ( [3TT i, we can obtain the following relation 



<T23 (0) 



1 . J32 - r23 . , (cos^(e/2)r23-sin^(e/2)r32) -2(r,3+r,,)/ 
- sin h sin 9 e ^ " '^' 

2 r23 + r32 r23 + r32 



~ sin0e-(™-^'*n'+^'-^+r^^)' (e'^' cos^ (6/2) - e-'" sin^ (0/2)) 



<T22 (0)) 



1 . ^r32-r23 . ^Sin2(0/2)r23-cos2(0/2)r32) -2(r„+r«)r 
- sm0 + sin0 e ■"'■'23+i32;f 

2 r23 + r32 r23 + r32 

+ ]- sin(6l)e"*^'=°'*' "n.+rjj+Fjj)/ ^^^ ^^^2 (g,/2) - g-''" sin^ (6»/2)) <T33 (0)> 

(sin"* (61/2) e-"=' + cos^ (0/2) e'^Oe"^™'' sn.+r^j+rjj)^ + 1 gjjj2 0g-2(r23+r32), 
+ ^ sin^ 0(e-2(r,3+r3.), _ ^-(cos= m,+r,3+r3,)r ^^^^^j^^ ^^^^ (q^^ ^ 



<T23 (0)) 



(32) 



Now, we obtain the evolution relation of the density matrix 
elements in the bare state representation. Since the expres- 
sions are very complex, here we only show the matrix ele- 
ments which will be used in the following. Based on these 
evolutionary matrix elements, we can write out the density 
matrix of the system in bare state representation at time t once 
the initial state is given, and then we can obtain the concur- 
rence of the density matrix. In what follows, we will discuss 
the entanglement dynamics and steady-state entanglement. 



A. Entanglement dynamics 

In the process of single-excitation energy transfer from the 
donor to the acceptor, the single excitation energy is initially 
possessed by the donor and the acceptor is in its ground state. 
Therefore the initial state of the system is 



Wm = \eg) = |772), 



(33) 



which means the initial conditions are that all matrix elements 
are zero except (t22(0)) = 1. According to Eq. ( l32l i. we know 
that the density matrix p{t) of the system belongs to the so- 
called X-class state. Then the concurrence can be obtained 
with Eq. 



at) 



1 

- sin 

2 



-32 



-23 



H-sin0 



r23 + r32 

(cos2(0/2)r 



23 



sin2(0/2)r32) 



-2(r23+r32)r 



— sinfle^^" 

2 



r23 + r32 

is-flni-i-r23+r32)/ 



: [e'" cos^ {912) - e''" sin^ (0/2))] 



(34) 



In what follows, we consider three special cases of interest: 
(1) The resonant case, i.e., wi = W2 = w„,, that is Aw - 0. 
Then the mixing angle 9 - ;r/2 and the energy separation 
e = 2^, thus we obtain 

1 



Cres(f) 



m2o 



(1 



.g-^(2flr') + ,-sin(^f)g4A'(2f)r' 



,(35) 



o 



! 0.5 



. 


■'-'"- 


h ;; ;;; ; ' 

", , , ; ' . '/ 


---r,„/7 = n.i - 

-■-T„./7 = 10 - 

— r,„/7 = ioo 


1 . .,>■ ' ;■ 


- 







7t 



FIG. 7: (Color online) The concurrence Crcs in Eq. i35t vs the scaled 
time yt for different bath temperature Ti„/y = 0.1 (dashed red line), 
10 (dash dotted blue line), and 100 (solid black line) in the resonant 
case Aoj/y = 0. Other parameters are set as y = 1, ^/y = 5, and 

AT/ J = 0. 



where N(2^) = ni(2^) + n2(2^) + 1. From Eq. ^, we find 
that the concurrence Cres(0 increases from zero to a steady 
state value 1 /N{2^) with the increase of the time t. Clearly, 
the steady state concurrence l/N{2^) decreases from one to 
zero as the temperature T„ increases from zero to infinite. In 
Fig.|7] we plot the concurrence (|35] | in the resonant case vs the 
scaled time yt and for different heat bath average temperatures 
Tm. Figure|2]shows the results as we analyze above. 

(2) The high temperature limit, i.e., T\,T2^> s. In this case, 
«i(e),n2(e) » 1, then we can have the approximate relations 
hi(e) =» «i(e) + 1 and h2is) ~ niis) + 1, which lead to r23 ~ 
r32. Then the concurrence ( |34] l becomes 



Chti(f) 



Sin(20) __,; „2 



sin^ 9N{e)yt _ ^^^ g^-{2 cos^ 8^+ \ sin^ 6N(E)y)t 



X [e""' cos^ (0/2) - e-'""' sin^ (0/2)) . (36) 

The expression of the concurrence (|36] | in the high tempera- 
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FIG. 8: (Color online) The concurrence Chti in Eq. l l36t vs the scaled 
evolution time yt for different mixing angle 9 = O.ln (dashed red 
line), 0.3n (dashed blue line), and 0.5n (solid black line) in the high 
temperature limit T„,/y = 100. Other parameters are set as 7 = 1, 
f/r = 5,^-1/7 =X2/r = O.Oir,,,, and Ar/y = 0. 



FIG. 9: (Color online) The concurrence Cm in Eq. J37t vs the scaled 
evolution time yt for different mixing angle 6 = O.ln (solid black 
line), 0.3;r (dash dotted blue line), and 0.5n (dashed red line) is plot- 
ted in the low temperature limit T^/y = 0.01. Other parameters are 
set as y = 1, ^/y = 5,Xi/y=X2/y = O.Oir,,,, and AT/y = 0. 



ture limit is not simple as that of the resonant case, but we can 
still observe the two points: The first is that the dependence 
of the concurrence on the angle is approximately sin 9; and 
the second is that the steady-state concurrence is zero, which 
means there is no quantum entanglement between the donor 
and the acceptor. This result can also be seen from the den- 
sity operator of the steady state for the donor and the acceptor 
In the high temperature limit, the steady state density matrix 
of the donor and the acceptor is p » {\eg}{eg\ + \ge}{ge\)/2, 
which is an unentangled state. Physically, this result is direct 
since the quantum systems will transit to classical systems in 
the high temperature limit. In Fig.[8l we plot the concurrence 
given by Eq. ( l36b vs the evolution time t for different mixing 
angles 0. Figure [8] shows that the concurrence experiences an 
increase from zero to a maximal value and then decreases to a 
steady state value with the scaled time -yt. 

(3) The low temperature limit, i.e., T ^ Q. Now we can 
approximately have h{s) ~ 0, which lead to r32 ~ and r23 ~ 
sin 6*7/2. Then the concurrence (|34] | becomes 

Citi(f) ~ sin 1 1 - 2 cos2(6»/2)e- ''""' ^^' + e^J '^' ^' 

X (e''^' cos^ (61/2) - e-'"" sin^ (6»/2))| , (37) 

where the subscript "Itl" stands for low temperature limit. 
Similar to the high temperature limit, the increase of the con- 
currence is also not simply exponential. The concurrence in- 
creases from zero to a steady state value sin Q with the increase 
of the scaled time f, which means the concurrence at long time 
limit is irrespective of the sign of the detuning. This long lived 
entanglement is much larger than that of the high temperature 
limit. We can also see the steady state concurrence from the 
viewpoint of quantum noise. When T x; 0, the steady state of 
the donor and the acceptor is p ^ |/l3)(/l3| with concurrence 
sin0. In Fig. in we plot the concurrence given by Eq. ((37) vs 
the evolution time t and the mixing angle 6. Figure |9] shows 



that the concurrence increases from zero to a steady state value 
with the scaled time t. 



B. Steady state entanglement 

From Eq. ( |34i ), it is straightforward to obtain the steady 
state concurrence between the donor and the acceptor. 



C. 



sinfl 



(38) 



In the high temperature limit, we have Chti(oo) ~ 0, and in 
the low temperature limit, we have Citi(oo) ^ sin 6*. For a 
general state, it is interesting to point out that the steady state 
concurrence Css depends on the temperature r,„ and the an- 
gle Q independently. For a given 6, the dependence on T^ 
is inverse proportional to N(e), and for a given r,„, the de- 
pendence on Q is sin0. In Fig. (TO] we plot the concurrence 
given by Eq. ( |38] l vs the temperature T^ for different mixing 
angles Q. Figure [10] shows that the steady state concurrence 
decreases with the increase of the temperature T,„. In Fig.[TT] 
we plot the concurrence given by Eq. ( |38T l vs the mixing an- 
gle Q for different average bath temperature T,„. Figure (TTj 
shows that the dependence of the concurrence on the mixing 
angle Q decreases with the increase of the average bath tem- 
perature r,„. Moreover, from Eq. ( [38] ), we can also see that 
the steady-state concurrence is independent of Ar at the high 
temperature limit. 



VI. CONCLUDING WITH REMARKS 

In conclusion, we have studied analytically coherent single- 
excitation energy transfer in a dimer consisting of a donor and 
an acceptor modeled by two TLSs, which are immersed in two 
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FIG. 10: (Color online) The steady state concurrence C.,., vs the bath 
temperature T,„ for different mixing angle = 0. l;r (solid black line), 
0.3 (dash dotted blue line), and 0.5;r (dashed red line). Other param- 
eters are set as y = 1, f /y = 5, and AT/y = 0. 
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acceptor decreases with the increasing of the bath tempera- 
ture. And the dependence of the steady state concurrence on 
the energy detuning is proportional to the sine function of the 
mixing angle and irrespective of the bath temperatures. More- 
over, we have found that the dependence of the steady state 
concurrence on the bath temperature difference is negligibly 
small with the current parameters. 

Finally, we give two remarks on the above obtained results: 
First, we should distinguish the present work from dynamic 
disentanglement suddenly or asymptotically (e.g., Refs. pS- 
dOII ). Mainly, there are three points of difference between 
the two cases: the initial state, the coupling between the two 
TLSs, and the coupling form between the TLSs and their heat 
baths. In dynamic disentanglement, the two TLSs is initially 
prepared in an entanglement state, there is no coupling be- 
tween the two TLSs, and the coupling form of the TLSs with 
their heat baths is off-diagonal. But in the present work, ini- 
tially the two TLSs are unentangled, there is a dipole-dipole 
interaction between the two TLSs, and the coupling form of 
the TLSs with their heat baths is diagonal. Certainly, the re- 
sults also differ In entanglement sudden death, the two TLSs 
disentangle to zero suddenly. But in this work, steady state 
entanglement is created. 

Second, in this work, we only address the problem about 
how is the dynamics of the created quantum entanglement in 
the process of excitation energy transfer |51]. But we do not 
address the question about the relation between initially pre- 
pared quantum entanglement among the pigments and the ef- 
ficiency for single-excitation energy transfer. Just as in quan- 
tum information science, quantum entanglement is considered 
an important resource since it can be used to enhance the ef- 
ficiency of quantum information protocols. Therefore it re- 
mains a question whether initially prepared quantum entan- 
glement can enhance the efficiency of excitation energy trans- 
fer 



FIG. 1 1 : (Color online) The steady state concurrence C,j vs the mix- 
ing angle 6 for different bath temperature T,„ly = 0.1 (dashed red 
line), 10 (dash dotted blue line), and 100 (solid black line). Other 
parameters are set as y = 1, ^/y = 5, and Ar/y = 0. 



independent heat baths. Special attention is paid to the effect 
on the single-excitation energy transfer probability of the en- 
ergy detuning and the heat bath temperatures of the two TLSs. 
It has been found that, the probability for single-excitation en- 
ergy transfer largely depends on the energy detuning in the 
low temperature limit. Concretely, the positive and negative 
energy detunings can increase and decrease the probability, 
respectively. In the high temperature limit, however, the effect 
of the energy detuning on the probability is negligibly small. 
We have also found that the probability is negligibly depen- 
dence on the bath temperature difference in the low and high 
temperature limits. We have also studied analytically quantum 
entanglement in the dimer system through calculating quan- 
tum concurrence. It was found that quantum entanglement 
can be created during the process of excitation energy trans- 
fer The steady state entanglement between the donor and the 
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Appendix: Derivation of quantum master equation Q 

In this appendix, we present a detailed derivation of quan- 
tum master equation (|7]i. Let us start from the Hamiltonian ([T]l 
of the total system. In the interaction picture with respect to 
the Hamiltonian 



Hq - //tlSs + Hb, 

the interaction Hamiltonian (|6]l becomes 



(A.l) 



Hiit) = 0-iiBii(f) + 0-22^22(0 + cr33B33(0 

+cr23e"'B23(0 + (Ti2e-'"Bl^{t), (A.2) 



11 



through introducing the following noise operators: 
A(t) + B{t), 



Bn{t) 

Biiit) 
Bait) = 



cos^{ei2)A{t) + sin^(0/2)B(O, 
sin2(6»/2)A(f) + cos^ {6 1 2) B{t), 



sin0 



\ k J 



,(A.3) 



with 



A{t) = e'"" 'AiQ)e-'"« 



= J]8ij(ay'^"'' + aje-"^'''), 



B(t) = e'"'^''B(0)e-'"'^'' 



YjS2k(b^ 



k^ 



+ bke-"^'"'). 



(A.4) 



Obviously, A(t) and B(t) are Hermitian operators. Note that in 
Eq. ( IA.2I ) we have made rotating wave approximation. 

Under the Born-Markov approximation, the master equa- 
tion reads 05l] 



Ps 



y-»CX 

Jo 



dTTrB[H,(t), [H,(t - T),ps ®Pb]1 (A.5) 



where Tr^ stands for tracing over the degrees of freedom of 
the heat baths. The density matrix pg = pjl^ ®P,/, of the heat 
baths means the two independent heat baths being in thermal 
equilibrium. 



(a) 
P,h 



where we denote the partition functions Za 



(A.6) 



7(«)^ 



7W^ 



TrBjexp(-/3ii/7)] and Zb = TrBjexp(-y62//^ ')] with 
ySi = l/Ti and jSa = ^ITi being respectively the inverse 
temperatures of the heat baths of the TLS 1 and TLS2. 

By using Eqs. ( IA.2l i and ( IA.5b and making rotating wave 
approximation, we can obtain the following quantum master 
equation 

PS = cr,2Pscr2i e''" {B23(-T)Bl^mdT 

Jo 

/-»oo 

+cr23Pscr32 e'"{Bl^i-T)B23mdT 

Jo 

/-*oo 

-T22PS e'"{B2mB\^(-T))dT 

Jo 



-0-33PS 



f 

Jo 



e-"'''{B\J0)B23{-T))dT 



'23^ 



+ 2_i ^minPsCTn 
m,n=l,2,3 



f 

Jo 



{B„„{-T)B,n,„mdT 



/-*oo 



«=1,2,3 

+h.c.. 



(A.7) 



where the correlation functions for the bath operators are de- 
fined as {X{t)Y{t')) = i:rB[X{t)Y(t')pB]. Notice that here we 
use the property (^(O^f')) = {X{t-t')Y(<d)) = {X(0)Y(t' - 1)) 
of the correlation functions for the bath operators. To de- 
rive the quantum master equation we need to calculate the 
Fourier transform of the correlation functions in Eq. (IA.7I) . 
For simplicity, here we only keep the real parts of the Fourier 
transforms of the correlation functions and neglect their imag- 
inary parts since the imaginary parts only contribute to the 
Lamb shifts, which are neglected in this work. According to 
Eqs. ( IA.3I ). ( IA.4b . and ( IA.6b . we can express the Fourier trans- 
forms for the correlation functions in Eq. ( IA.7b as follows: 



f 

Jo 







dT{Bn(0)Bn(-T)} 



dT{B22iO)B22(-T)) 



dT{B33iO)B33(-T)) 



f 

Jo 



dT{Bu(-T)Bn(0)r 



Da + Db, 

dT{B22(-T)B22(0)r 

Jo 

cos^{e/2)DA + ^m\ei2)DB, 



f 

Jo 



dT{B33(-T)B33(0)r 



/-*oo 

dT{B22(-T)Bum 

Jo 

f 

Jo 



sin^{e/2)DA + cos\ei2)DB, 



-f 

Jo 



dT{Bu(-T)B22m 



= cos2(6'/2)D* + sm^(0/2)Dl, 

/-*CO 

dT{B33(-T)Bum = dT{Bll(-T)B33m 

Jo 
= sin^(e/2)D'^ + cos^(0/2)Ds, 

dT{B22(-T)B33m 

Jo 

ym\0/2) (Dl+Dl), (A.8) 



dT{B33i-T)B22m = 



where the parameters are introduced as 



Da = dT{A(0)A(-T)), Db^ dT{B(0)B(-T)), 

Jo Jo 



Dl - f 

Jo 



dT{A{-T)A{0)} 



^ Dl^ f 

Jo 



dT{B{-T)B(Q)}. 

(A.9) 



Since the noise operators B„„ (n - 1, 2, 3) are Hermitian oper- 
ators, then we have the relations 



f 

Jo 



f 

Jo 



dT{B„„{-T)B„„m ^ dT{B,„,n(0)B„„(-T)} 



(A. 10) 



Therefore we can know all the diagonal correlation functions 
in Eq. ( |A.7b as long as we obtain the expression of Da and 
Db- According to Eqs. ( IA.4b and ( IA.6I ), we can calculate the 
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expression of Da as follows, 



1 

1 



= Re[DA] + /Im[D^], 



(A. 11) 



where 



Rq[Da\ = lim ;r£)iMgiM[2ni(w)+ 1]. (A.12) 



Note that in the third line of Eq. dA.llI ) we have used the for- 
mula: 



f 

Jo 



dre"^ 



■ n6{u)) + iP—, 

0) 



(A. 13) 



where the sign "P" stands for the usual principal part integral. 
Similarly, we can obtain the expression of Re[Z)s], 

Re[£)s] = Um ngi (w) gj (w) [2«2 (w) + 1] . (A. 14) 

<jj-»0+ 



Here gi(w) and £)2(w) are the densities of state of the heat 
baths of the TLSl and TLS2, respectively. And ni(w) = 
l/[exp(J5ia)) - 1] and h2(cj) = l/[exp(/?2w) - 1] are the av- 
erage thermal excitation numbers. Using the same method we 
can obtain the following expressions: 



Re 



Re 



and 



dTe'"{B23(0)B\^{-T)) 
Jo 

/-*oo 

dTe-''\B23(-T)Blm 
Jo 



n 
4 



= ■- sin2 6[gi(e)g](e)(hi(e) + 1) + e2(e)gl(e)(«2(e) + 1)] 



(A. 15) 



Re 



Re 



/-*oo 

dTe-''^Bl^(0)B23i-T)) 
Jo 



f 

Jo 



dTe"^{BU-T)B23m 



= - sin^ e[gi{e)gj(s)ni{s) + g2{s)gl(s)n2(s)UA.l6) 



By substituting these correlation functions into Eq. (IA.7l l and 
returning to the Schrodinger picture, we can obtain quantum 
master equation O. 
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